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ABSTRACT 

Context. Studies of dust evolution and processing in different phases of the interstellar medium is key to understanding the lifecycle 
of dust in space. Recent results have challenged the capabilities and validity of current dust models, indicating that the properties of 
interstellar dust evolve as it transits between different phases of the interstellar medium. 

Aims. We characterize the dust content from the IC 434 H ii region, and present a scenario that results in the large-scale structure of 
the region seen to date. 

Methods. We conduct a multi-wavelength study of the dust emission from the ionized gas, and combine this with modeling, from large 
scales that provide insight into the history of the IC 434/Ll630 region, to small scales that allowe us to infer quantitative properties 
of the dust content inside the H ii region. 

Results. The dust enters the H ii region through momentum transfer with a champagne flow of ionized gas, set up by a chance en¬ 
counter between the L1630 molecular cloud and the star cluster of cr Ori. We observe two clearly separated dust populations inside 
the ionized gas, that show different observational properties, as well as contrasting optical properties. Population A is colder (~ 
25 K) than predicted by widely-used dust models, its temperature is insensitive to an increase of the impinging radiation field, is 
momentum-coupled to the gas, and efficiently absorbs radiation pressure to form a dust wave at 1.0 pc ahead of cr Ori AB. Population 
B is characterized by a constant [20/30] flux ratio throughout the Hii region, heats up to ~ 75 K close to the star, and is less efficient 
in absorbing radiation pressure, forming a dust wave at 0.1 pc from the star. 

Conclusions. The dust inside IC 434 is bimodal. The characteristics of population A are remarkable and can not be explained by cur¬ 
rent dust models. We argue that large porous grains or fluffy aggregates are potential candidates to explain much of the observational 
characteristics. Population B are grains that match the classical description of spherical, compact dust. The inferred optical properties 
are consistent with either very small grains, or large grains in thermal equilibrium with the radiation held. Our results confirm recent 
work that stress the importance of variations in the dust properties between different regions of the ISM. 

Key words. Keywords should be given 


1. Introduction 


The advent of high-resolution infrared (IR) imaging during the 
last decade has provided indisputable evidence that dust resides 
within the ionized gas of Hii regions, rather than being associ- 
ated with photo-dissociation regions alon g the line of sight (e.g., 
[Povich et al.||2Q07t [Paladini et al.||2Q12| ). More than often, the 
discussion in literature is focussed on the question of what kind 
of dust population is able to survive inside an H ii region, as the 
physical conditions in these regions predict that dust will be re¬ 
moved and/or destroyed in the ionized gas, either through radia- 
tion pressure, a stellar wind, dust sublimation, or dust sp uttering 
( |Inoue|2002t [Everett & Churchwell|2Q10[|Draine|201l| ). 

The work of Chini et al. ( 1986| ) triggered the thought of 
a separate dust component inside Hii regions, which was sup¬ 
ported by the model presented in [Everett & Churchwell ( |2010[ ), 
who discussed the possibility of evaporating cloudlets overrun 
by the ionization front (IE), resupplying the ionized gas with a 
new generation of dust. [Paladini et al.[ ( |2012[ ) showed that the IR 
emission from evolved H ii regions is typically s tratified as fol- 
lows: polycyclic aromatic hydrocarbons (PAHs; [Tielens 2008[ ) 
and cold dust (~ 20 K) trace the PDR and a dense shell of the H ii 
region, respectively, while 24 pm emission peaks inside the ion¬ 
ized gas, close to the central source. Ochsendorf et al. ( 2014b[ ) 
proposed a new scenario to simultaneously explain the presence 
and morphologies of the mid-IR emission inside Hii regions. 


Whereas several authors attribute the 24 pm emission inside ion¬ 
ized regions to the increase of VSGs compared with BGs ( [Elagey[ 
[et al.|2Qllt [Paradis et al.|2011[ ), the analysis of the IR arc in IC 
434 (Ochsendorf et al.|2014^ from now: Paper 1) revealed that 
the increased heating by stellar photons from the nearby star can 
explain the mid-IR emission. 


Besides the intrinsic nature of the dust inside Hii regions, 
recent studies ( Salgado et al. ( [2012 ), Paper 1) have indicated 
that dust properties inside H ii regions may differ from that used 
in widely-used dust models of the diffuse interstellar medium 
(ISM) ( |Li & Draine 200 ![ [Compiegne et al.[[2011[ ). In addi¬ 
tion, results from the Planck mission indicate that dust properties 
evolve when transiting from one region of the ISM to another 
(e.g., Planck Collaboration et al.|2014a[ ). Here, we focus on dust 
grains entrained within photo-evaporation flows, and follow the 
evolution as they transit from the molecular cloud phase to ion¬ 
ized regions of space. 


We examine the dust content and processes governing the 
properties of the dust inside H ii regions, using the IC 434 region 
as a benchmark model. The multitude of observations available 
for the region, its large angular extent, and its profitable location 
well outside the Galactic plane, makes IC 434 a suitable candi¬ 
date to pursue this question. We extend the work and analysis 
described in Paper 1 to the far-IR. The new observations allow 
us to probe the dust content of the IC 434 region on a wide range 
of wavelengths which reveals a new dust population in the H ii 
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region. When combined with the theoretical model described in 
Paper 1, the analysis provides a means to constrain the proper¬ 
ties of the dust entrained in the ionized gas flow of the IC 434 
region. In particular, we will argue that the dust content in IC 
434 is best explained through a bimodal distribution, which dif¬ 
fers signiflcantly from that seen in the diffuse ISM. We discuss 
our findings in the light of dust evolution inside the molecular 
cloud from which the dust is evaporated, such as the formation 
of aggregates through coagulation, and through recent findings 
of dust evolution as constrained by, e.g., the Planck mission. 

In Sec.[^ we set the stage by digging into the history of the 
large-scale structure of the IC 434/Orion B molecular cloud re¬ 
gion, that provides us with an updated and improved model for 
the ionized gas flow in IC 434 (as compared to Paper 1), neces¬ 
sary to derive the properties of the dust introduced into the H ii 
region. In Sec. we present the observations. After character¬ 
izing the global properties of the dust, we combine the obser¬ 
vations with detailed modeling in Sec.[^ to derive properties of 
the emitting grains coming from the ionized gas. The results of 
the modeling are described in Sec.[^ We discuss our findings in 
Sec|3 and conclude in Sec. [ 8 ] 


2. Observations 


The observations used in this work are in part described in Paper 
1, but are extended with Herschel observations at T > 70 /im, 
Planck observations ( [Planck Collaboration e t al.|2014b ), and th e 
Planck dust model R1.2Q ( [Planck Collaboration et al.|[MT4^ . 
PACS 160 yum and SPIRE 250 /im, 350 jam and 500 //m pho¬ 
tometry were taken from the Herschel Science Archive (HSA) 
from the Gould Belt survey ( Andre et al.||2010 ) with observa¬ 
tion IDs 1342215984 and 1342215985. Herschel data was inter- 
calibrated with IRAS 70 jam (PACS 70 //m), Planck 857 GHz 
(PACS 160 /im, SPIRE 250 /im and 350 /im), and Pl anck 545 
GHz ( SPIRE 500 /im) through the method described in |Bemard| 
jet al.|(|2010|). The Herschel resolutions are 8.4” (70 /im), 13.5” 
(160 /im), 18.2” (250 /im), 24.9” (350 /im), and 36.3” (500 /im). 


respectively, for parallel mode observations at (60” s“^) scan 
rates. 

The reduction of Spitzer/IRS observations is described in 
Paper 1, but is extended here with the Short-Low (SL) modules, 
covering the 5 - 15 /im region. The IRS SL detector suffered 
from a background level, causing apparent gradients throughout 
the modules and mismatches between the SLl and SL2 spectra. 
This background level persisted after background subtraction of 
the BCD images and affected the reduced spectra. Eollowing 
Sandstrom et al.j ( |2012| ), we opted to correct the BCD images 
with the use of the inter-order regions that resulted in a signifi¬ 
cant improvement of the quality of the spectra. 


3. The expanding HI shell surrounding Orion OB1b 

In this section, we connect the IC 434 H ii region with the large- 
scale structure of the Orion Belt region. We will argue that the 
large IR bubble, approximately centered on the open cluster of 
(T Ori, is blown by the Belt star e Ori, and that cr Ori has crossed 
the distance from the Orion OBlc region, currently approaching 
the LI630 molecular cloud. Besides that the analysis provides a 
new insight into the history and evolution of the Belt region, the 
argumentation ultimately leads to the development of the flow 
model described in Sec. and, subsequently, the derivation of 
grain properties that are entrained within the ionized gas. 


3.1. The GS206-17+13 shell: mass and kinetic energy 


The HI shell GS206-17-rl3 ( jEhlerova & Palous|2005| ), centered 
at (l,b) = (206.1,-17.5), has a central radial velocity of vlsr = 
12.8 km s“^ and an expansion velocity of roughly Vgxp = 8 km 
s“^ It extends ~ 2° x 4° (13 x 26 pc at 385 pc; see below) 
in HI, and is also seen at mid- to far-IR wavelengths (Eig. [^. 
Apparently, it is centered on the Orion OB lb association and cr 
Ori. The total mass in the shell is estimated by adopting a dust 
emissivity per unit mass of gas, Ky, which is often modeled as 
a power law in the form of Ky = (lO)v^, with the frequency v 
in THz ( Beckwith et al.|p^90| and p the spectral index of the 
dust opacity. It should be noted that this relation could be more 
complex, depending on the exact grain properties (e.g., [Planck 
[Collaboration et al.] ( [2014a[ ), and references therein). We use the 
t?> 5 ?> (optical depth at 353 GHz) and the p maps from the Planck 
dust model R1.20 ( [Planck Collaboration et aL||2014a[ ), and de¬ 
fine an aperture that encloses the bubble shell seen in WISE and 
Planck images. We estimate the contribution of Galactic emis¬ 
sion along the line-of-sight unrelated with the bubble shell by 
determining the value of T 353 in a region close by at the same 
Galactic latitude (~ 23% of the total emission), and subtract this 
from the optical depth measured within the aperture. With Ky = 
2.1 cm“^ g“^ at 353 GHz, and yS = 1.5 as constrained by the 
Planck model, we calculate the total dust mass of the shell 
= (^ 353 /^y) 5 ', with S the physical size of the emitting region. We 
And a total dust mass of = 34 M© or, equivalently, a gas mass 
of Mg = 3400 Mq (using a gas-to-dust ratio R = 100). By dis¬ 
tributing the mass back into a homogeneous sphere with radius 
20 pc (the average radius of the HI shell), we And an average 
density of ^locai = 4.1 cm“^, i.e., the original density in which the 
bubble expanded. Then, the total kinetic energy of th e shell can 
be estimated with the expansion model described in [Chevalier 
( [T?74l ): 
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Here, r^h is the observed radius of the shell. Using ^locai =4.1 
20 pc and Vexp = 8.0 km s“\ we And E’tot = 5 x 10"^^ 


cm 


^sh ■ 


erg, which is only a factor of few lo wer than the typical M netic 
energy injected by a SN (~ 10^^ erg; Veilleux et al. ( |2005 )). 


3.2. Driving source of the GS206-17+13 shell 

The GS206-17-rl3 shell is filled with ionized gas that is known 
as the IC 434 emission nebula. The bright gas emission results 
from the ionizing flux of the central massive component of the 
cr Ori cluster, the triple-star-system cr O ri AB at J = 385 pc 
(Simon-Dfaz et al.[2OTT[[Caballero[2008| that is part of a quintu¬ 


ple system ( [van Loon & Qliveira|2b03| ' Caballero|2007a ), evap 
orating the LI630 cloud and launches a champagne flow into the 
expanding bubble. In this respect, it is tempting to attribute the 
expanding shell to the continuous input of the stellar wind of cr 
Ori AB to its surroundings. However, cr Ori is a young system 
(4ge = 2-3 Myr; Caballero][2008[ ) which e xhibits a weak-wind 
(log(M) = -9.7 M© yr“^ Vcx, = 1500 km s~^; Najarro et al.| 2 ()TT ). 


This amounts to an integrated mechanical luminosity of Lmech 
= 1/2 M 4ge = 1.3 X 10^^ erg, of which only ~ 10% will 


couple as k inetic energy to the ISM ( [van Buren[[1986[ [Veilleux 
et al.[2005| ). This excludes cr Ori as the driving source of the 
expanding GS206-17-rl3 shell. 

The Orion OB lb sub-association is cent ered on Orion’s belt 
and has age estimates ranging from 1.7 Myr ( Brown et al.|1994 ) 
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Fig. 1. The WISE mid-IR view of the region surrounding the Orion OB association at 12 jjm. The different sub-groups of the Orion 
OB association are overplotted for stars of spectral type B2 and earlier ( [Brown et al.|1994| ), where the oversized symbols represent cr 
Ori (diamond) and e Ori (inverted triangle), respectively. The 12 yum emission reveal the edges of the Orion molecular clouds as they 
are being illuminated by the ionizing flux of members of the Orion OB association. Also shown are several prominent structures in 
the image: the LI630 molecular cloud, the emission nebula NGC2024, and the Orion Nebular cloud (ONC). Of particular interest 
here i s the bubble centered at (/,/?) = (206.1,-17.5), which is the IR counterpart of the GS206-17-rl3 HI shell ( jEhlerova & Palous 


2005| ) and encircles the members of Orion OB lb (dashed outline). Overplotted in solid lines are the space motions of the stars. 


tracing back the motion of the stars in the plane of the sky during 1 Myr. Eor cr Ori, this trajectory has been extended to 3 Myr, 
where the open diamond symbols mark 1 Myr, 2 Myr, and 3 Myr. 


to 8 Myr ( |Blaauw|1964| ). However, as noted by |Bally| ( |2008| ), the 
supergiant status of the Orion Belt stars (f Ori, e Ori, and 6 Ori) 
imply an age of the lb association of at least 5 Myr. We find 11 
stars of spectral type B2 or earlier within the Orion OB lb asso¬ 
ciation ( [Brown et al.|1994| ). Remarkably, the massive members 
of the Orion OB lb association are positioned along the GS206- 
17-rl3 shell (Eig.[^, which provide clues to their origin. 



the kinetic energy of the GS206-17-I-13 shell within ~ 11 Myr 
(assuming a constant wind strength and 10% kinetic efficiency), 
which is a reasonable timescale compared to the estimated age 
of the Orion OB lb association. Surely, the supergiant stars f Ori 
and 6 Ori could have been part of the formation of the GS206- 
17-rl3 shell. However, the Hipparcos dist ances towards C Or i 
and 6 Ori, d = 225^28 ^ “ ^^^-23 P^ (van Leeuwen|2007), 


respectively, compared to the average distance of the stars lo¬ 
cated within the velocity extent of the expanding HI shell (vlsr 
= 12.8+7.8 km s~^), d = 440 pc, renders this scenario implau¬ 
sible. This difference in distance is reflected in the low radial 
velocity vlsr = 1 km s“^ for both f Ori and S Ori, compared 
to the velocity extend of the HI shell, which places the stars out¬ 
side of the bubble. Perhaps, the stars have been thrown out of the 
clu ster through dynamical inte raction in the young OB lb clus¬ 
ter. Reynolds & Ogden] ( [1979[ ) already noted that 6 Ori does not 
have a clear H ii region surrounding it, implying that the ionizing 
radiati on can escape freely into the Orion-Eridanus superbubble 
cavity ( Brown et al.|1995 ) instead of illuminating the surround¬ 
ings of the star, i.e., the GS206-17+13 shell. We conclude that 
the stars at d = 440 pc are the best candidates for driving the 
expansion of the GS206-17+13 shell, with e Ori as the main 
contributor. 
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3.3. Space velocities of Orion OB1b and approach ofcr Ori 4. Dust emission from the IC 434 region 


We calculate LSR space velocities of the stars belonging t o the 
Orion lb association following the method described in |Cox| 
et al.| ( |2012| ), using proper m otions and radial velocit ies, and 


correcting for the solar motion ([Cogkunoglu et al.|201 1 ). Proper 
motions are taken from |van Leeuwen| ( |2007| ), except for the case 
of cr Ori, where due to the complex nature of the system the new 
reduction of the Hipparcos catalogue did not produce an accept¬ 
able solution {pa cos^ = 2 2.63+10.83 mas and ps = 13.45+5.09 
mas, [van Leeuwen||2007| versus p^cosS = 4 .61+0.88 mas and 
Ps = -0.40+0.53 mas, Perryman et al. 1997| ). The new reduc¬ 
tion is also inconsistent with proper motion measurements from 
other members of the central quintuple system of cr Ori, cr Ori 
D (paCosS = 0 .79+0.92 mas and ps = -0.90+0.47 mas, |van| 
Leeuwen||2007| ) and cr Ori E (pa cos S = 2.2+0.92 mas and ps 


= -1.4+0.47 mas. 


Caballero|2007b|). He nce, for cr Ori the proper 


motion values of |Perryman k al. { 1991) are used. Figure[2over- 


plots the resulting space velocities for stars with vlsr within the 
velocity extent of the expanding GS206-17+13 HI shell. The 
plotted trajectories can be traced hack in time, and show that it 
is possible that the stars have dispersed to their current location 
projected along the bubble shell. In particular, the stars at the 
northern side of the GS206-17+13 shell at longitude I > 206° 
seem to be moving out of the bubble, influencing the morphol¬ 
ogy of the shell traced at 12 pm in the direction of movement. 
This asymmetry is also seen as an elongation of the HI shell in 
Galactic longitude ( [Ehlerova & Palous|2005| ), and indicates that 
the Orion OB lb stars indeed have a relative velocity with respect 
to the Orion B molecular cloud and the GS206-17+13 shell. We 
add here as a caveat that due to the the location of Orion in the 
anti-direction of the Galactic center, the motion of the OB as- 
sociation is mostly directed radially away from the Sun ( |Brown| 
et al.||1994[ |de Zeeuw et al.| 1999 1, which complicates precise 


ation (> 5 Myr, Sec. |3.2[ ). Therefore, [Bally| ( |2008| ) tentatively 
classifled the cr Ori cluster as being part of the Orion OBlc as¬ 
sociation. Following Paper 1, the radial veloc ity vt sr of cr Ori 


and the LI630 cloud are similar (~ 10 km s ; Maddalena et al 


19861 [Gibb et al.||1995| ). The proper motion of a Ori translates 


to a velocity of 10 km s“ , we will use this value as the rela¬ 
tive velocity , Vreb between the cloud and cr Ori AB in Sec. 
Then, the trajectory of space motion at 1 Myr, 2 Myr, and 3 Myr 
shows that cr Ori AB can cross the distance from the Orion OB Ic 
region to its current location. Indeed, the sub-ass ociation Orion 
OBlc has an estimated age of 2 M yr to 6 Myr ( [Blaauw| 1964[ 
[Warren & Hesser|1978[|Bally|2008| ), which would solve the age 
discrepancy as described above. However, this assumes that the 
observed space velocity represents the relative velocity between 
cloud and star, but the actual traversed distance could be smaller 
if the bulk motion of the Orion region is co-moving with the 
stars (i.e., Vrei < 10 km s“^). In any case, as the cr Ori cluster ap¬ 
proaches the LI630 cloud, the ionizing flux reaching the cloud 
surface increases, resu lting in the bright photo-evaporation flow 
seen to date (Sec. |5.2| ), known as the IC 434 emission nebula. 


measurement of proper motions. In addition, the distance to¬ 
wards some of the more distant stars in Orion OB lb gives rise to 
Hipparc os parallax errors of 40 - 50%. In this re spect, the Gaia 
mission ( jde Bruijne|2012 Perryman et al.|2()0T ) is expected to 
improve proper motion measurements and distance determina¬ 
tions towards the members of Orion OB 1. 

Based on its spatial proximity, cr Ori has often been classi¬ 
fled as being part of the Orion OB lb association. jBallyj (|2Q08|) 
alread y noted that the age of the system (~ 2-3 Myr; jCaballeroj 
|2008| ) is incompatible with the age of the Orion OB lb associ- 


4.1. Morphological structures 

Figure ^ shows the mid-IR view of the region surrounding IC 
434. The WISE 12 pm is dominated by PAH emission tracing 
the PDR surfaces directly adjacent to the IF. However, the gas 
content in the same region shows a different morphology (Fig. 
1^): the bulk of the Ha emission concentrates within a smaller 
region with radius ~ 0.75° (~ 5 pc at a distance of 385 pc) inside 
the mid-IR bubble (Fig.|^). The ionized gas originates from the 
evaporation of the LI630 molecular cloud by the ionizing flux 
of cr Ori AB. This gas flow effectively introduces dust into the 
ionized gas through momentum transfer. 

Figure shows the distribution of dust inside the IC 434 re¬ 
gion at wavelengths 8 pm > A> 160 pm. The 24 pm dust emis¬ 
sion surrounding cr Ori AB has earlier been attributed to a dust 
wave (Paper 1). Her e, we report the detection of emission from 
the ionized gas (Sec. |4.3| ) at 8 pm, 12 pm, and at Herschel wave¬ 
lengths longwards of 160 pm. The 12 pm WISE Alter is broad 
(ranging from 7 pm -17 pva), and in principle is sensitive to both 
warm dust radiating at 100 K - 200 K and transiently heated par¬ 
ticles, such as Very Small Grains (VSGs) and IR fluorescence 
from PAHs. While PAHs mainly show up in broad emission fea¬ 
tures at 5 pm -18 pm, VSGs typically reveal a continuum from A 
>13 pm, which could be a real co ntinuum or a superposition of 
multiple blended spectral features (jCesarsky et al.|2000l|Abergel 
jet al.|2002j ). 

After inspection of the images at different wavelengths, one 
can distinguish two separate morphological IR emission struc¬ 
tures inside the Ha emitting volume (Fig. [3). On the one hand, 
there is a component seen at all wavelengths but for 24 pm (mor¬ 
phology ’A’). This component is characterized through its large 
diffuse scale and clearly associated with the full extent of the 
LI630 cloud. The morphology is spread throughout the entire 
IC 434 nebula, ultimately culminating in an arc-shaped increase 
in emission to the east of cr Ori AB at a projected distance of 
J ~ 1 pc (best observed at 8 pm), together with a clear drop¬ 
off in emission to the west of the star. For the remainder of this 
work, we will refer to this arc-like structure at J = 1 pc as the 
outer dust wave (ODW), due to the striking resemblance of the 
structure with the dust wave observed at 0.1 pc (see below). On 
the other hand, there is the emission detected with WISE and 
Spitzer at 22 pm and 24 pm which is clearly distinct from mor¬ 
phology A, peaking at J ~ 0.1 pc and revealing a relative narrow 
flow (morphology ’B’). This flow seems to be connected to a 
bright part of the LI630 cloud which is almost symmetrically 
positioned around the Horsehead Nebula, again forming an arc¬ 
like structure that peaks in intensity at J ~ 0.1 pc. We will refer 
to this structure as the inner dust wave (IDW). Close to the star, 
the IDW lights up at 12 pm and 8 pm due to heating of the grains 
and an increased dust-to-gas ratio (Paper 1). 


4.2. Spectral energy distributions 

To get insight in the properties of the emitting dust in the H ii 
region, we constructed SEDs from the WISE, MIPS, PACS and 
SPIRE photometric points. Ah images were convolved to the 
point spread function of the SPIRE 500 pm band using the ker¬ 
nels described in j Aniano et al. j ( |MTT] ) . Subsequently, a common 
’dark’ spot in the maps was chosen for background subtraction 
(centered at 05^37"^00^ -01^52"^290. The absolute uncertain¬ 
ties on the photometry are estimated at 10% for the WISE and 
Spitzer images, while we use uncertainties of 20% for PACS and 
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Fig. 2. (A) Two color composite in the mid-IR of the region around cr Ori AB (white asterisk). Green is WISE 12 jam, whereas red is 
WISE 22 jjm. The large IR counterpart of the GS206-17+13 HI shell surrounds cr Ori AB. The east-side of the bubble is marked by 
the large molecular cloud L1630, which includes the Horsehead Nebula (HH). The projected distance between the Horsehead and 
cr Ori, as shown by the w hite line along which the Spitzer/IRS observations are positioned (Eig.[^, is about 0.5 degrees or 3.7 pc (at 
385 pc; |Caballero 2008). The white cross reveals the offset position used for the Spitzer/IRS background subtraction. (B) Smoothed 
Ha image from SHASSA (4’ resolution). The gas emission around cr Ori AB (IC 434) originates from a photo-evaporation flow of 
the LI630 molecular cloud. The bright part appears roughly circular with a radius ~ 0.75°. (C) Schematic view of the same region, 
showing the IR bubble (solid outline) and the Ha emission (dashed circle), created by the photo-evaporation flow (arrows) of the 
LI630 molecular cloud (filled black rectangle). 


Morphology A 


Morphology A-hB 


Morphology B 




Fig. 3. Two different infrared morphologies within IC 434. In all panels: north is up, east is to the left. Morphology A is evident at 
every wavelength except for 24 jjm, and shows a diffuse structure emanating from the entire LI630 molecular cloud. Morphology 
B is a relative narrow structure, ranging from the LI630 molecular cloud towards cr Ori AB (white asterisk in the 160 jam images), 
where it interacts with radiation pressure from the star at projected distance 0.1 < J < 0.4 pc (Paper 1), forming the inner dust wave 
that is visible at wavelengths T < 24 yum. The 8 yum and 12 yum images show a combination of both morphologies (A-rB). The lower 
panels zoom into the direct surroundings of cr Ori AB and cover the region in the dotted white box seen in the upper panels. The 
8 yum, 12 yum and 160 yum images show another arc-like structure at a distance of ~ 1 pc away from the star (dashed line), forming 
the outer dust wave. The 12 yum suffer from several ghost features (’G’). The 160 yum image reveals an extended collimated flow 
(dotted line) away from the cr Ori AB, not seen at 8 yum and 12 yum, which is possibly an evaporating member (’E’) of the large cr 
Ori cluster. The color bars show the brightness scale in units of MJy/sr. 
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Fig. 4. (a) The mid-IR view of the Horsehead - cr Ori AB (yellow asterisk) region. Red is MIPS 24 jjm, blue is Ha from the KPNO 
4.0m Mayall telescope. Overplotted in white rectangles are the regions for which Spitzer IRS/spectra and SEDs were extracted. The 
gray rectangles are the location of the peak intensities of the outer dust wave (ODW) and inner dust wave (IDW) for which only 
the SED is available, (b) Erom left to right, Spitzer/IRS SL spectra at a distance J = 3.7 pc (located partly inside the Horsehead 
PDR), 3.1 pc, 2.0 pc, 0.8 pc (inside the Hii region) and 0.3 pc (projected on top of IDW) from cr Ori AB are extracted, where d is 
the distance along a straight line connecting cr Ori AB and the Horsehead. Indicated are gas emission lines, the PAH plateau at 6 - 
9 jam and at 11 - 14 jam, and the incident radiation field Gq. PAH emission is clearly observed throughout the Hii region, (c) Same 
as (b), but for the IRS LL modules. Labelled beneath the curve is the source of the rising continuum in the mid-IR. (d) Ear-infrared 
spectral energy distribution of the regions. In all cases except for region 7, a single-component fit is conducted on the SED, as VSGs 
emit significantly at T < 70 jam. In region 7, we assume thermal equilibrium for the warm component and add a second component 
to the fit. (e) Spectral energy distributions for the ODW and IDW. 


15% for SPIRE, respectively ( [Bernard et al.||2010| ). We extract 
SEDs corresponding to the regions covered by IRS (regions 1- 
7;Fig.|g and complement these with two regions at both dust 
arcs, located at 1 pc (ODW) and 0.1 pc (IDW), respectively. The 
results are shown in Eig.|^ 


4.3. Are we tracing emission from the ionized! gas? 

The analysis of the remainder of this paper assumes that the dif¬ 
fuse dust emission from morphology A and B (Sec. |4.1| ) origi¬ 
nates from within the ionized gas of the IC 434 emission nebula. 
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For morphology B, this is obvious from the bright interaction 
zone close to cr Ori AB (Paper 1). Below, we will argue that the 
same is true for morphology A. 

Both for the continuum and the spectroscopic observations, 
we have subtracted a background value to exclude dust asso¬ 
ciated with the GS206-17-rl3 shell and contamination along 
the line-of-sight (Fig. [^. This method may leave some resid¬ 
ual structure if the bubble shell is not homogeneous. However, 
H 2 emission at 9.7 yum, 12.3 yum, and 17.0 yum immediately dis- 
sapear (Fig. |^&c) when moving from region 1 (dominated by 
the Horsehead PDR) towards region 2-7, indicating a transi¬ 
tion from the shielded PDR to the atomic/ionized gas, attest¬ 
ing to the success of the correction for the background emission 
from any surrounding neutral material. Moreover, morphology 
A is clearly separated from the LI630 molecular cloud, hhing 
the Hof emitting region (Fig. with diffuse emission, ranging 
from the LI630 molecular cloud to cr Ori. Furthermore, mor¬ 
phology A contains a limb-brightened arc at a projected distance 
of 1 pc east of the star (more pronounced in associated 8 yum 
emission; see below), upstream in the flow of ionized gas, and 
it abruptly terminates when moving towards the west-side of the 
star, both indicative of an interaction with cr Ori AB, and argu¬ 
ing against a connection to background material. Indeed, in the 
photo-evaporation flow, most of the material is expected to be lo¬ 
cated near the symmetry axis connecting the star and the cloud 
(Sec. |5.2| ). Within this framework, the total diffuse dust emission 
located between the star and the cloud represents a dust reservoir 
that has evaporated of the LI630 molecular cloud, and is moving 
along with the ionized gas, some of which has not yet reached 
the interaction zone ahead of cr Ori AB. 

Emission observed in the WISE 12 yum band and IRAC 8 
yum closely follow the spatial morphology of the cold dust of 
Morphology A. As can be seen from the Spitzer/IRS spectra 
shown in Fig.|^&c, the WISE 12 yum is picking up PAH emis¬ 
sion from the 7.7 yum and 11.3 yum band. The dust continuum 
from VSGs inside the Hii region (region 2-6) only rises at X 
> 20 yum (see Sec. |4.3.1| ), well outside the WISE hlter edge (~ 
17 yum). We note that emission from [Nell], arising from the IC 
434 emission nebula, potentially contaminates the detection in 
the broad WISE hlter. However, if the [Nell] emission were to 
dominate the observed hux in the WISE band, this would fol¬ 
low the smooth, exponential appearance from the Ha emission 
measure (Fig. |^), which is not seen in the observations. We 
conclude that WISE traces PAH emission. The IRAC 8 um band 
also traces PAH emission. By realizing that the cold dust ob¬ 
served by Herschel, and the PAH emission as seen by WISE and 
Spitzer are spatially related (Fig.|^, we can use the evolution of 
the mid-IR PAH spectrum to answer the question if morphology 
A traces a dust component that is located inside the ionized gas 
(Fig.§. In this respect, we note that jCompiegne et aL] ( |2007| ) al¬ 
ready showed that the PAH emission originates from within the 
ionized gas through correlation with ionized gas tracers. 

The intensity ratio of the 7.7 yum to 11.3 pm PAH fea¬ 
tures is a well-known indicator of the charge state of PAHs 
(e.g., Szczepanski & Vala|1993 Bauschhcher|2002| ), as neutral 
species emit signihcantly less in the 6 - 9 pm region compared 
to charged species ( [Ahamandola et al.|1999| ). Figure plots the 
1 .ll\\3 pm intensity ratio, which initially decreases from region 
1 to 2: |Compiegne et al.| ( |2007| ) argue that, initially, the PAHs be¬ 
come more neutral due to the increase of free electrons across the 
IF. Subsequently, from region 2 to region 6 the 7.7/11.3 pm ratio 
shows a steady increase inside the H ii region, indicating an over¬ 
all increase of the (cationic) charge state of the PAHs. Thus, the 
PAHs are located inside the ionized gas: we are probing emission 


10^ 



Fig. 5. The evolution of the PAH 7.7 yum over 11.3 yum intensity 
ratio as a function of distance into the H n region, tracing the 
ionization state of the PAH molecules (e.g. |Bauschhcher||2002[ 
Szczepanski & Vala|1993|). The intensit ies of the bands are mea- 
sured with PAHFIT ( [Smith et al.|2007| ). The ratio shows a sys¬ 
tematic rise with distance into the H ii region, indicating that the 
molecules reside within the ionized gas, and a dramatic drop in 
the region closest to the ionizing star. The numbers correspond to 
the regions drawn in Fig.[^ the dotted lines show the Icr uncer¬ 
tainty constrained by the quality of the spectra . The upper axis 
plots the ionization parameter y ( [Tielens [2005 1 ). 


from regions that are located increasingly closer to the star. The 
combination of a decrease of electrons available downstream in 
the ionized flow wit h a gain of stel lar photons increases the ion¬ 
ization parameter y ( [Tielens|2005 ), deflned as jn^. Here, 

Go is the hardness of the radiation held in units of the Habing 
held ( Habing[[l968 ), T is the temperature of the (ionized) gas, 
and Hq is the electron number density. A ’typical-sized’ inter¬ 
stellar PAH with 50 carbon atoms will become more positively 
charge d as y increase s, and the cationic state will dominate at y 
>10"^ ( Tielens|20Q5| . In region 7, we are tracing a region PAH- 
free cavity around the star and we cannot trace the charge ra¬ 
tio reliably anymore (see Sec .[7.3[ ). The observed relation of the 
7.7/11.3 pm ratio with y (Fig.[^ is in reasonable agreement with 
observations of the Orion Bar and NGC 2023 ([Galliano et al.[ 
[2008] ). 

Another proxy for the charge state of PAHs is the 11.0 pm 
featu re due to C-H out of plane bending modes of cationic PAHs 
(e.g., [Hony et al.|20'0H [Rosenberg et al.|2011 ). Indeed, we And 
a hint of the 11.0 pm feature in regions 5-6 close to the star, 
but the quality of the IRS spectrum prevents us from making 
a conclusive statement on a detection of the 11.0 yum satellite 
feature. 

In summary, we argue that morphology A represents a dust 
population entrained within ionized gas of the IC 434 emission 
nebula because of the following: 

(1) We subtracted the background to exclude dust associated 
with the GS206-17-rl3 shell: analysis of the Spitzer/IRS 
spectrum reveals the absence of the H 2 feature immediately 
after leaving the Horsehead PDR: there is no neutral region 
associated with the spectra of regions 2-7 (Fig.|^. 

(2) PAH emission shows a tight correlation with the long wave¬ 
length dust emission associated with morphology A (Fig.[^. 

(3) PAHs are increasingly charged when moving to regions pro¬ 
jected closer towards the cr Ori AB. Thus, the PAHs are 
located ins ide the ionized gas. Thi s conclusion was also 
reached by [Compiegne et (2007), based upon observa¬ 
tions of the Horsehead nebula. Consequently, the same holds 
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Region 

1 

2 

3 

4 

5 

6 

7 

Location 

PDR 

Hii 

Hii 

Hii 

Hii 

Hii 

DW 

log(Go) 

2.0 

2.2 

2.3 

2.5 

2.9 

3.3 

4.2 

[20/30] 

0.28 

0.20 

0.21 

0.19 

0.22 

0.21 

0.33 


of VSGs in the H ii region, which contribute to the emission ob¬ 
served at /I < 70 yum, this is not a straightforward task (Anderson 
et al.|2012] ). As the SED clearly reveals the presence of a warm 


and the [20/30] continuum flux ratio for the regions in Fig. 
1^ Region 1 is located partly inside the Horsehead photo¬ 
dissociation region (PDR), whereas regions 2-7 are located in¬ 
side the H II region. Region 7 is projected on top of the dust wave 
surrounding cr Ori AB. The [20/30] ratio is insensitive to the in¬ 
crease of Go inside the H ii region, except in region 7, where the 
dust reaches thermal equilibrium (see text). 


for the cold dust emission from morphology A that is seen 
co-spatial with the PAH emission. 

(4) The diffuse emission from morphology A culminates in an 
arc-shaped emission around a Ori AB, accompanied by a 
dear-drop in emission and a cavity behind the star, both in¬ 
dicating an interaction of the dust with cr Ori AB close to the 
symmetry axis connecting the star and cloud (Sec.[^. 

4.3.1 . VSGs: 15 yum - 35 yum mid-infrared spectra 

The IRS observations were designed to follow the evolution of 
the mid-IR spectrum in an increasing radiation field. The 15 yum 
- 35 yum spectra (Fig. [J) show gas emission lines from [S III] 
and [Si II], typical for an ionized nebula. In addition, as region 
1 is (partly) located inside the Horsehead PDR, molecular H 2 
emission is seen along with the PAH emission plateau, ranging 
from 16 yum -19 yum. A steep continuum at T > 15 yum reveals the 
presence of VSGs, as the intensity of the radiation field, log(Go) 
= 2.0, is too low for classical grains, or Big Grains (BGs), to 
produce continuum emission in the mid-IR. This continuum de¬ 
creases between region 1 and region 2, but persists throughout 
the HII region and shows a remarkable constant spectral slope 
(see Tab. [T]) as revealed by the flux continuum ratio at 20 yum 
and 30 yum, [20/30]. The carbonaceous nature of the carriers of 
the mid-IR continuum emission (often assumed for the VSGs 
in dust models) is confirmed through the absence of the silicate 
emission feature at 2 0 yum (e.g.,|van Boekel et al.|2003 ), but for 
an opposing view see Li & Draine ( |2002| )~ 

The constant slope strongly implies the presence of 
stochastically-heated VSGs that survive throughout the Hii re¬ 
gion, as a constant flux ratio does not reflect a dust population in 
TE inside an increasing radiation field. The decrease in observed 
flux between the PDR (region 1) and the Hii region (region 2 - 
7) likely comprises a sub-population of VSGs that do not sur¬ 
vive the transition from the PDR into the H ii region. We have 
labelled the responsible carriers in Fig. [^: note that in region 
7, where log(Go) = 4.0, both VSGs and BGs are expected to be 
in TE and can both be responsible for the increase of the ratio 
[20/30] in region 7. In this respect, the observed dust tempera¬ 
ture of this component at the IDW (73-82 K), and the reduced 
UV-to-IR opacity (Sec. 4.3.4[ ) argue in favor of BGs, consistent 
with the findings in Paper 1. Note that these BGs can both be of 
carbonaceous or silicate nature; the peak temperature is too low 
for the 10 yum and 20 yum silicate features to be in emission. 


dust component (Fig.|^, we exclude the 24 yum and 70 yum data 
for all regions, except for region 7 and the IDW. Here, we assume 
that the dust component emitting at 24 yum an 70 yum has reached 
TE with the radiation field, and perform a two-component black- 
body fit to the SED. 

The cold component reveals a striking constant appearance 
throughout the H ii region, characterized by temperatures of 
~ 23 K - 26 K, and low spectral emissivity, yS ~ 0.7 - 1.0. Stiil, 
these values of p are within what is expected for the derived 
dust temper atures ( [Veneziani et al.||2010t Planck Collaboration 
et al.||2011| ). Note however, that the parameters and p seem 


degenerate, and derived values are sensitive to systematic effects 
induced from photometric calibration uncertainties to different 
grain temperatures superpositioned along the line-of-sight (see 
|Juvela & Ysard|20T^ , and dust optical properties ( |Jones|20l^ . 
Nevertheless, the constant appearance of the SED is curious, 
considering that the impinging radiation field varies significantly 
(Tab.[^. For the warm component, we fix y(5 = 1 due to the lack 
of data points with respect to free param eters, a value appropri - 
ate for the spectral index in the (mid-)IR ( [Mennella et al.|19^ ). 
Then, the warm component reaches a maximum temperature of 
7'd,H = 82 K at the position of the IDW (J = 0.1 pc). 

As the cold dust does not contribute much at 70 yum at the 
location of the IDW (< 10%), the addition of the extra compo¬ 
nent at long wavelengths does not change the results obtained 
in Paper 1. However, the peak dust temperature of the hot dust 
depends on the choice ofyS: for a value ofyS = 1.8 for the hot com¬ 
ponent (as was used in Paper 1), the peak temperature decreases 
to rd,H = 73 K. The addition of the hot component does affect the 
fit in region 7 by increasing p and lowering because here 
the hot component contributes significantly at 70 yum. This im¬ 
plies that the single-component fits in region 1-6 overestimate 
T^ c slightly. Indeed, by fixing the s pectral emissivity dXp = 1.5 
(as mea sured by Planck for IC 434 Planck Collaboration et al 
(2014a|)), yield = 19 K - 20 K for region 1-6. The larger 
value of p as measured by Planck is most likely caused by differ¬ 
ent dust populations and mixing of temperatures, as the Planck 
model samples the spectral emissivity at a 30’ resolution, and 
does not correct for contamination along the line-of-sight. The 
high angular resolution of Herschel, and our careful background 
subtraction isolates the emission from the ionized gas. Still, to 
determine the actual temperature and spectral emissivity of dust 
in the H ii region with a single MBB fit is not straightforward, 
due to the presence of multiple dust components inside the ion¬ 
ized gas, as reflected by the SEDs. 

In Sec. \4A\ we showed that the far-IR emission in IC 434 
is dominated by morphology A, which fills the H ii region and 
forms the ODW at J = 1 pc to the east of the star. The dust popu¬ 
lation accountable for morphology A is distributed inside the H ii 
region, ranging from region I (d = 3.1 pc) to the ODW (d = 1 pc) 
(Sec. |4.3| ), close to the symmetry axis (Sec.|^. Radiation pres¬ 
ume from the star prohibits the dust from drawing closer (Sec. 
|5.3| ). In a first attempt to constrain the properties of the grains 
responsible for morphology A, we use DUSTEM ( [Compiegne 
et al.|20lf] ) to predict the shape of the dust FIR SED by using a 


4.3.2. Cold dust: far-infrared spectral energy distribution 

It is common to model the observed dust SED with a single mod¬ 
ified blackbody (MBB) function. However, due to the presence 


radiation field evaluated at the projected distance to cr Ori AB. 
Compared to the single-temperature MBB fit, DUSTEM allows 
for a range of equilibrium temperatures, a natural consequence 
of the size distribution and composition of interstellar dust. 

Figure shows the predicted SED by DUSTEM for both re¬ 
gion 1 and the ODW. The result are shown for two scenarios. 
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Reg. 1 

ODW 

IDW 


rf(pc) 

log(Go) 

3.7 

2.0 

1 

3.0 

0.1 

5.1 

MBB 

Ti,c (K) 

23 (19) 

27 (20) 

82 (74) 

DUSTEM 

req,s (K) 0.1 lim 

T'eq.s (K) 1.0 jum 

Teq.c (K) 0.1 jim 
T'eq.c (K) 1.0 jum 

31 

20 

34 

21 

47 

31 

55 

33 

100 

65 

130 

74 


Table 2. Dust temperatures for grains at the location of re¬ 
gion 1, the outer dust wave (ODW), and the inner dust wave 
(IDW), derived through 1) a modified-blackbody-fit (MBB) and 
2) DUSTEM modeling of the equilibrium temperature for a 0.1 
yum silicate grain and 1.0 yum silicate grain (Teq^s) and amorphous 
carbon grain (Teq^c)- Also listed are the projected distan ce J, and 
the hardness of the radiation field, Go ( |Habing|[l968| ). The ob¬ 
served cold component T^ c is calculated with yS as a free param¬ 
eter (see Fi g.[4|), and by choosing yg = 1.5 as gi ven by the Planck 
dust model ( [Planck Collaboration et al.|2014a| ), denoted between 
brackets. 


First, we use grain species (both amorphous carbon and silicates) 
and size distributions that successfully reproduce the extinction 
and emission properties of a selected region of the diffuse ISM 
at high Galactic latitude (the DHGF model; jCompiegne et al.j 
[201 Ij ). In this case, the small end of the size distribution domi¬ 
nates in surface area, which causes the predicted dust emission 
to peak at short wavelengths, causing DUSTEM to be unable 
to reproduce the observed SED. Second, by only including very 
large silicate grains in the size distribution a reasonable fit can 
be obtained, as explained below. Table compares T^ c from 
the modified-blackbody fit, at the location of region 1 and the 
ODW, with equilibrium dust temperatures of graphite and sili¬ 
cate grains calculated with DUSTEM. The results for the IDW 
region are also listed. While DUSTEM predicts that the grains 
are significantly heated, resulting in a total temperature increase 
of ~ 50% between region 1 and the ODW (because of the in¬ 
crease of stellar fiux as they approach the star), the modified 
black body fits reveal that the cold dust component only in¬ 
creases ~ 15% in temperature. This result highlights the main 
puzzle presented by the FIR SEDs of the Hii region: a cold 
dust population, whose temperature is rather insensitive to an 
increase in radiation field. Below, we will outline two scenar¬ 
ios that might explain the remarkable properties of the grains 
responsible for Morphology A. 

A reasonable fit for region 1 can be obtained by only includ¬ 
ing silicate grains larger than 1 micron to the fit. Silicates have 
lower equilibrium temperatures compared to amorphous carbon, 
and the size distribution composed of very large grains pushes 
the peak of the predicted dust emission to longer wavelengths. 
Still, this grain population is not able to fit the observed SED at 
the ODW region, as the predicted dust emission peaks blueward 
of the observed SED. This can perhaps be solved if the grains de¬ 
viate from the classical spherical description. For example, the 
equilibrium temperature Tgq of fluffy aggregates is thought to be 
lower, comp ared to their compact counterparts (typically about 
10% - 20%; |Fogel & Leung|1998| ). This is because UV, visible 
and near-IR absorptivity remain similar as fiuffiness is increased, 
whereas the sub-mm emissivity increases significantly ( jBazell &] 
Dwek|1990] ). An intuitive explanation for the presence of large 


aggregates inside the H ii region arises when realizing that the 
grains are ’freshly’ evaporated from the molecular cloud (Sec.[^ 
in which the grain are thought to coagulate to large, fluffy aggre- 



Wavelength (/i-m) Wavelength (/xnn) 

Fig. 6. Modeling of the far-infared spectral energy distribution 
using DUSTEM ( jCompiegne et al.|201 1| ). Two panels are shown, 
corresponding to two different regions (Fig. [^. The left panel 
shows region 1 (J = 3.7 pc); the right panel displays the outer 
dust wave region (ODW; d = I pc). Overplotted is the predicted 
dust emission for two different composition and size distribution 
of grains (see text). 


gate structures ( jOssenkopf & Henning|1994 Ormel et al.|2009|). 
This argumentation connects well with the work of Martin et al. 
( 2012| ), that used BLAST and IRAS images of selected regions 
in Vela of moderate column density to show considerable varia¬ 
tions in, e.g., the submillimeter opacity (as compared to the dif¬ 
fuse ISM). The authors argued that the inferred grain properties 
may reflect past histories in the evolution of the grains, such as 
coagulation in dense regions. 

Besides a rather rigorous change in the size distribution and 
invoking fiuffiness or porosity, another explanation for the ob¬ 
served low dust temperatures for Morphology A (as compared 
to predictions made by DUSTEM) can be linked to recent re- 
sults from, e.g., the Planck mission. Planck Collaboration et al. 
( |2014a| ) argue that Tobs is not a simple tracer of the radiation 
field, but rather probes variations in dust properties (grain struc¬ 
ture, material change, and size distributions). Similar results 
are reported in [Plank Collaboration et al.| ( |2014| ): through com¬ 
parison with the jPraine & Li ~ pQQ7[ ) physical model, the au¬ 
thors show that the location of the peak wavelength of the dust 
SED does not only trace the intensity of the radiation field (like 
DUSTEM), but also intrinsic variat ions of dust properties. Other 
results from the Planck mission ( |Planck Collaboration et al.| 
2014a[ ) include that dust in specific region of the ISM is more 
emi ssive than those used in du st models from jLi & Draine[ ( |2001[ ) 
and jCompiegne et al.j ( 2007] ), which would lower the equilib¬ 
rium temperature with respect to that comput ed using these grain 
models. In conclusion and as pointed out by Jones ( 2014j ), dust 
models that do not allow for the evolution of grain properties as 
they transit between different regions of the ISM, will be chal¬ 
lenged in the light of new observational evidence obtained by 
Planck and Herschel. 


4.3.3. Dust- and gas mass inside IC 434 

The total dust mass in IC 434 is determin ed by deriving the av- 
erage T 353 from the Planck model R1.20 ( [Planck Collaboration 
et al.|2014a| inside an aperture that encloses the bright Ha emis¬ 
sion from the emission nebula. We find T 353 = 1.5 x 10“^. The 
total dust mass can the n be estimated directly through M^jr = 
—S (jHildebrandl 19 8 3 1) , where S is the size of the aperture and 

If 

Kx is the grain opacity at wa velength A. If /C 353 = 1.9 cm^ g“^ 
( jWeingartner & Draine|2001j ), then the total dust mass in IC 434 
equals M^jr = 2.3 Mq. 

The total gas mass inside IC 434 can be estimated by in¬ 
tegrating the gas density profile, given by the Ha observations 
(Sec j5.2j), out to a distance of 9 pc (the apparent extent of the 
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bright Ha emission). This yields a total hydrogen column of A^h 
= 1.6 X 10^^ cm“^. We assume that the gas evaporates of a 9 pc 
X 9 pc surface from the LI630 molecular cloud (Sec. HI- Then, 
the total gas mass in the IC 434 emission nebula is Mg 100 
Mq. The total dust-to-gas ratio in the IC 434 region is therefore 
~ 0.02, similar to that seen in the diffuse ISM (~ 0.01). A thor¬ 
ough study of gas-to-dust ratios in large-scale photo-evaporation 


nected these observations to the molecular cloud phase from 
which the grain are evaporated. However, even these large grains 
cannot explain the temperature dependence that seem insensi¬ 
tive to an increase in intensity of the radiation field. This might 
confirm recent results, indicating that the peak wavelength of 
the SED is not determined by the intensity of the radiation field 


flows could reveal if dust entrainment significantly enriches the 

et al. 

[2014 

1, but also traces changes in grain properties (e.g.. 

ISM with grain materials. 

Jones|2014 



4.3.4. Dust opacities of the ODW in the UV and IR 

Following the method and using parameters described in Paper 
1 for the IDW, we estimate the optical depth of the ODW at UV 
wavelengths through aperture photometry at Tuv = 1.2 x 10“^, 
with a corresponding dust mass of Md,uv = 7.4 x 10“^ Mq. In 
cont rast, fo llowing the same method for Mh tu as described in 
Sec.|4.3 3 and by using a:i6o = 10 cm^ g“^ |Weingartner & Drainej 
|2001| ), then the total mass in the ODW from the IR observations 
at 160 yum is M^jr = 0.16 Mq. Thus, we find for the dust content 
in the ODW, M^jr/M^^uv = 22. This effect was already noted for 


the warm co mponent in Paper 1 and in the W3 region ( [Salgado 
|et al.pof^ , and implies that the dust opacities at UV and IR 
wavelengths for the cold dust component are also very different 
from that seen for the diffuse ISM: this may well be a general 
characteristic for dust inside Hii regions. Similar results are ob¬ 
served in (molecular) clouds where the dust opacity is altered, 
reflected in a flattening of the dust extinction curve towards UV 
wavelengths ( [Cardelli et al. 19891. Often, this is attributed to 
dust coagulation or a dearth of small particles in the dust size 
distribution. 


4.3.5. Summary on the observations 

The total dust-to-gas ratio inside IC 434 is ~ 0.02. The ODW 
contains 5 - 10% of the total dust content (2.3 Mq) inside IC 
434. The dust content inside the ionized gas can be divided in 
two different populations, following two separate morphologies 

(Fig.g. 

Morphology B (Fig.|^ shows the signatures of carbonaceous 
VSGs in regions 1-6, responsible for the constant [20/30] flux 
inside the H ii region. The emission from this component is most 
pronounced at 24 jjm. The reduced UV-to-NIR opacity, and the 
inferred peak dust temperature (73 K - 82 K) at J = 0.1 pc indi¬ 
cate that large grains also populate morphology B and perhaps 
dominate the emission from the IDW (results from Paper 1). 
These BGs are only visible close to the star as the grains heats 
up. 

Morphology A (Fig. is a combination of cold dust at T > 
160 yum, and mid-IR emission from PAHs at 8 yum and 12 yum 
emission that is tightly correlated with the cold dust throughout 
the Hii region. The PAHs show an increasing level of charge 
state in regions at smaller projected distances from the ionizing 
star, providing convincing evidence that the PAHs are located 
inside the ionized gas and therefore, the grains from morphology 
A are located in the H ii region as well. 

Dust temperatures of morphology A are between 23 K < 

< 26 K (depending on the exact value of JS); this is colder than ex¬ 
pected for a typical size distribution of dust in the diffuse ISM, as 
predicted by DUSTEM. We have pointed out that the cold tem¬ 
perature can be explained by (porous/fiuffy) micron-sized grains. 
This is consistent with the reduced UV-to-IR opacity compared 
to that observed for dust in the diffuse ISM, and we have con- 


5. The flow of dust into the IC 434 region 

As has been discussed in Sec. |3.3| cr Ori approaches the L1630 
molecular cloud to create a champagne flow that enters into the 
expanding bubble, blown by the Orion OB lb association (Fig. 
[^. The increase in the ionizing flux impinging the LI630 molec¬ 
ular cloud with time will strengthen the champagne flow. Note 
that this scenario differs from Paper 1, that considered a station¬ 
ary champagne flow. The ionized gas accelerates into the H ii re¬ 
gion and, subsequently, drags the dust along with it. The work in 
Paper 1 describes the theoretical framework to derive the proper¬ 
ties of the dust grains introduced into the H ii region and provides 
an analysis of the IDW surrounding cr Ori AB. In this section, 
we extend this work by considering the emission from morphol¬ 
ogy A as well (observed at 8 yum, 12 yum, and long wards of 160 
yum). Furthermore, we present how one can use dust waves to 
constrain geometry (porosity/fiuffiness) and the radiation pres¬ 
sure efficiency of the grains. 


5.1. Modeling the flow of ionized gas 


Henney et al.| ( |2005} investigated the evaporation of a cloud illu¬ 
minated by a point source, and the structure of the ionized flow 
emanating of the cloud surface. In their hydrodynamical models, 
the photo-evaporation flows are found to be steady during a large 
part of their evolution, in which the flow properties change only 
on timescales much longer than the dynamic time for flow away 
from the IF. In the steady phase, the flow resembles a ’cham¬ 
pagne flow’ ( |Tenorio-Tagle|| 197 9| ) or a ’globule flow’ ( Bertoldi 
1989| ), depending on the curvature of the IF. 

An IF that engulfs a fiat, homogenous cloud illuminated by 
a point source will become concave, and the photo-evaporation 
flow emanating from this geometry will resemble a champagne 
flow structure, as th e streamlines do not diverge effectively 
( [Henney et al.||2Q05] ). In practice, the ionized gas from a con¬ 
cave IF will accelerate more slowly compared to , e.g., the photo - 
evaporation flow from a convex-shaped globule ( |Bertoldi|1989| ). 
The relative slow acceleration of the gas will increase the effec¬ 
tive flow thickness, /zgfF, defined as n^hQ^ = f dr, where no 
is the ionized gas density and r is the distance from the star to 
the IF. Thus, the shape of the IF (convex, fiat, or concave) and 
the incident ionizing flux uniquely determine the structure of the 
resulting photo-evaporation flow. 

We envision the approach of cr Ori to the LI630 molecular 
cloud as follows. We consider the LI630 molecular cloud as a 
fiat IF, ionized by an (approaching) point source. By combin¬ 
ing the density prescription that reproduces the observed emis¬ 
sion from the IC 434 region (Eq. wit h the results of the nu¬ 
merical models of Henney et al.| ( [2005] ), we developed a flow 
model that improves the model described in Paper 1 substan¬ 
tially. We restrict the modeling of the flow to the symmetry 
axis connecting the star and the cloud, because the flow prop¬ 
erties are dominated by the immediate surroundings of the sym- 
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met ry axis ( Henney et al. |2005 ), as will be addr essed in Sec. 
|5.2[ M oreover, the photo-evaporation models from |Henney et al.| 
( |2005| ) show that in the case of a flat IF illuminated by a sin¬ 
gle point source, the streamlines do not diverge at the distance 
between star and cloud, so we can consider the gas flow to 
be moving one-dimensionally (i.e., orthogonal to the surface of 
the molecular cloud) constant along the line-of-sight. Thus, the 
symmetry axis provides a decent approximation for the (two- 
dimensional) structure of the flow between the cloud and the 
star. In short, the analysis results in an increased line-of-sight 
extent of the cloud compared to Paper 1, but is consisted with 
the numerical models of |Henney et al.| ( [2005] ) and comparable to 
the size of the LI630 molecular cloud. 

The photo-evaporation flow from the LI630 molecular cloud 
resembles a champa gne flow following the prescription of 
Tenorio-Tagle \\919) . In this scenario, the density contrast be¬ 


tween a cloud and the inter cloud medium drives a shock into the 
inter cloud medium, while a rarefaction wave travels back into 
the cloud, setting up a champagne flow that accelerates the gas 
to supersonic velocities as the gas moves away from the IF. The 
density distribution of the cham pagne flow can be described by 
an exponential density gradient (|Bedijn & Tenorio-Tagle|| 19811 
|Tielens|2QQ5] ): 


n{r) = no exp 


r-R. 


Ro 


( 2 ) 


J = Jo 


Jo 


Mn(r)fdr = ^^^(1 - exp [-2/?s//?o]), (3) 

2ct 


where Jq = Qol^Trr^ is the ionizing flux from the star, in the 
case that none of the photons are absorbed within the H ii region. 
Equation [^yields the solution: 


(4) 



Here, no is evaluated at the IF and equals J/c^, where J is the 
ionizing flux reaching the cloud and is the isothermal sound 
speed in the ionized gas. The value R^ is the distance between the 
cloud and the star, and Ro = c^t is the scale length of the flow, 
where t is the time since the ionized gas started to stream into the 
Hii region. We follow the conclusions of [Henney et al. ( |2Q05| ) 
by adopting a constant flow thickness as the star approaches the 
cloud, i.e., the scale height of the flow, Ro (see Eq.|^, will adjust 
itself to the distance between the IF and the star, R^. In p articular, 
the effective flow thickness for a flat IF is ^ 0.35/?s (jHenneyj 
jet al |2QQ5| Fig. 7). 

Thus, assuming that the density gradient can be described by 
Eq.|^ (Paper 1) and an effective flow thickness of /leff = 0.357?s, 
one can solve for the scale height Ro of the exponential flow 

by writing n^ Hqq = ' (n(r))^ dr. To determine the complete 

structure of the flow as the ionizing star draws closer, we now 
only have to solve for the ionizing flux, J, reaching the cloud 
surface at any given time. The ionizing flux will increase as the 
star advances towards the cloud and will accelerate the IF, in¬ 
creasing the ionized gas density streaming into the H ii region. 
Consequently, the number of recombinations between the star 
and the cloud will increase and, in its turn, lower the amount of 
photons reaching the cloud surface. The actual amount of pho¬ 
tons reaching the cloud is determined by integrating the number 
of recombinations in a column r extending from the star (r = 0) 
to the IF (r = R^): 


This result is similar to the one reported in Bedijn & Tenorio- 
Tagle| ( T98T] ), except their solution is written in terms of number 
density; we opted to follow the method from jSpitz^ ( j 197 8| ) and 
provide the solution in term of stellar fluxes. With the density 
prescription, we can solve for the scale height Ro of the flow 

using the previous mentioned relationship, n^ /zefF = Jq ' (^(^))^ 
dr. With no = //Cs, and = 0.357?s, it is easily shown that Ro 
= 0A6R,. 


5.2. Flow parameters 

Consider the LI630 molecular cloud as a uniformly distributed, 
extended cloud. We take for the ionizing flux of cr Ori AB 
log(2o) = 47.56 photons s“\ appropriate for an 09.5V star 
( [Martins et al.||2005] ). We let t he star approach the cloud at a 
velocity of 1 0 km s~^ (Sec. 3.3), which after 3 Myr (the adopted 
age of cr Ori; [Cabahero[2008| iresults in the conflguration seen to 
date: a projected distance of 3.7 pc between the LI630 cloud and 
cr Ori AB. At r = 0, we assume that a steady photo-evaporation 
flow already Alls the medium between the star and cloud (see 
Fig.[3].. 

Figure [7] reveals the evolution of the flow parameters, Vf and 
zzh, as cr Ori trav els towards the cloud at a velocity of v* = 10 
km s“^ (Sec. |3.3[ ). The velocity of the flow, Vf, is obtained after 
integration of Eq.[^ resulting in a linear velocity law. Th e initial 
veloc ity of the gas is Vf = 10 km s“^ (a D-critical front; Spitzer 
1978 ) and, because of the constant scale height Ro ( Henney et al 
2005), reaches the same maximum velocity of Vf = 32 km s“^ at 
ah times when the gas reaches the location of the star. Taking 
into account the relative velocity between star and cloud, this 
amounts to a maximum relative velocity of 42 km s“^ 

The density structure of the flow aX t = 3 Myr reaches n^ = 
35 cm“^ at the cloud surface. This density is offset by a factor 
of ~ 3 with the earlier derivation described in Paper 1. This dif¬ 
ference follows from the choice of Qo'. in Paper 1, the 09.5V 
model atmosphere of [Schaerer & de Koter[ ( [1997] ), with an ion- 


4nr'^2cl 


(l-cxp[-2RJRo]). 


izing luminosity of log(2o) = 48.25 photons s“\ was taken to 
calculate the incident flux on the LI630 cloud. Here, we use in¬ 
stead the updated model atmospheres from [Martins et al.[ ( [20Q5[ ) 
and take log(2o) = 47.56 photons s“^ The decrease of ionizing 
photons effectively lowers n^. In addition, we calculate the bal¬ 
ance between recombinations and photo-ionization appropriate 
for a champagne flow structur e with Eq. [4| wherea s in Paper 1 
we had followed the work from [Abergel et al. ( 2002[ ) and Spitzer 
( 1978[ ) that calculate Jo in a spherical geometry. While a spher¬ 
ical, globule-like geometry is appropriate for the small evapora - 
tion zone for the Horsehead considered in [Abergel et al. ( [2002[ ), 
it clearly is not applicable to the entire champagne-flow from the 
L1630 molecular cloud (e.g.. Fig . [7]). The S[III] lines that were 
used in [Compiegne et ah ( 2007[ ) (and that we had adopted in 
Paper 1) are close to or below the critical density, and its usage 
to derive reliable densities is questionable. Therefore, lacking 
other reliable density tracers that can pin down the initial den¬ 
sity at the IF, calculating the amount of photons that reach the 
molecular cloud in time, using the appropriate density laws from 
Eqs.[^and[^ is the most reliable way in constraining the flow 
parameters in the IC 434 nebula. In addition, we compare with 
numerical models from [Henney et al. ( [2005 ) to obtain a better 
description of the ionized flow of the LI630 molecular cloud. 

The lower density nu remains compatible with the observed 
Ha emission measure from the IC 434 region, as the electron 
dens ity dep ends on the scale size of the emitting gas through n^ 
= ^JEM|l. Here, EM and / are the emission measure and the 
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exponential density gradient inside the H ii region, depending on 
the choice of Rq. We look into this dependency by investigating 
our results after varying the scale height Rq in Sec. [7] 

In summary, we have updated the photo-evaporation model 
compared to the one described in Paper 1 by using (1) a more 
recent model atmosphere, (2) a better prescription of calculating 
the ’insulating’ layer of recombining material streaming of the 
molecular cloud (Eqs.[^&[^, (3) a comparison with numerical 
results for photo-evaporation flows from [Henney et al.| ( |2005| ), 
and (4) and allowed for the time evolution of the flow as cr Ori 
approaches the cloud. The exponential density gradient (Fig.[^) 
and the adopted scale height (Sec. reproduces the observed 
Ha emission measure well if the IC 434 region measures 9 pc 
along the line of sight, in agreement with the photo-evaporation 
models of |Henney et al.| ( [20n5] ). Note, that while the updated flow 
structure does affect some of the results derived in Paper 1, the 
main conclusions do not change and the derived grain properties 
are only affected marginally, as we will discuss in Sec.[^ 


5.3. The formation of dust waves 



Fig. 7. (a) Density structure n of the champagne flow as cr Ori 
moves toward the cloud (Sec. 3.3). The different curves represent 
the density structure after time t: 0 Myr, 1 Myr, 2Myr, and 3 
Myr. The dotted line shows the evolution of the density at the 
IF, uq. The grey dashed box shows the outline of the region of 
panel c. (b) Velocity structure of the flow. To transform this to the 
frame between star and cloud, one has to add the relative velocity 
between the star and cloud, = 10 km s“^ (c) The density 
profile 3Xt = 3 Myr, overplotted with the observed density profile 
from the Ha emission measure, assuming a constant scale size 
of the emitting gas along the line of sight of / = 9 pc (see text). 


scale size along the line-of-sight of the emitting volume, and 
(= ^h) is the electron density of the (fully) ionized gas. Thus, by 
lowering Hq by a factor of 3, / has to increase factor of 3^ along 
the line-of-sight to reproduce the EM of the ionized flow, which 
amounts to / = 9 pc at the LI630 molecular cloud, similar to the 
angular extent of the IC 434 region (Fig.p|). 

While strictly speaking, Eqs.[^and|^old for the symmetry 
axis only (extending from the star to the cloud radially), the bulk 
of the emission measure originates from the region within one 
unit distance R^ from the symmetry axis ( [Henney et al.||20Q5" 


Fig. 5a) and, therefore, the dust emission from the Hii region 
will largely be located within this region. In the case of IC 434, 
this amounts to an emitting area of ~ 7.5 pc centered on the 
symmetry axis, consistent with our adopted line-of-sight extent 
of the LI630 molecular cloud. Even though our adopted geom¬ 
etry and flow structure is consistent with numerical results from 
[Henney et al.| ( |2005] ) and the size of the LI630 cloud on the sky, 
we note that the observed EM can be used with any I to match an 


In this section, we briefly review the dust wave phenomenon. 
Following Paper 1, the equation of motion for dust is written as: 


mdVd 


JVd 

dr 




+ Icr^kT Hi 


8 ^ 

3 




1/2 


(5) 


The first term on the right hand side is the radiation pressure 
force, Frad, and contains: the geometrical cross-section of the 
grain, (t^\ the radiation pressure efficiency (averaged over the 
incident radiation field), and the speed of light, c. The sec¬ 
ond term is the collisional drag force ( [Draine & Salpeter|1979[ ), 
Fdrag, and contains: the Boltzmann constant, k\ the temperature 
and number density of the gas (of species /), T and ni\ and 
Si = with Vdrift the drift velocity of the dust 

through the gas. Furthermore, is the dust mass and Vd is the 
dust velocity. 

As the grains flow along with the gas into the Hii region 
through the drag force, Fdrag, the radiation pressure of the star, 
Frad, will act on a dust grain contained within the flow and will 
cause it to lose momentum. As the dust is slowed down, it will be 
pushed through the gas with drift velocity Vdrift, transferring mo¬ 
mentum towards the gas. We consider collisional drag only, as 
Coulomb interactions do not seem to dominate the drag term for 
dust-gas momentum transfer inside IC 434 (Paper 1). Eventually, 
for a dust grain approaching the star radially (i.e., zero impact 
parameter), the dust will be stopped at the equilibrium radius, 
Teq, upstream with respect to the star, where |Frad/Fdragl equals 
unity. If at this point the amount of momentum transfer between 
gas and dust is insufficient for the gas to be dynamically per¬ 
turbed, gas and dust will decouple significantly, resulting in an 
increase of dust upstream relative to the star as the gas flows 
along unhindered: the formation of a dust wave. Dust grains 
approaching the star with a non-zero impact parameter will be 
pushed around the star, resulting in an arc-shaped structure. For 
a detailed description of the physics of a dusty photo-evaporation 
flow we refer to (Paper 1). 


5.4. The equilibrium- and minimum radius of a dust wave 

The radiation pressure opacity, determines the trajectory of 
a dust grain as it moves toward and passes the star : 

= CtpCd/wd = QipC (6) 
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tQ = 0 Myr Iq ^ 2 Myr 



Fig. 8. The trajectories of two grains with different value of ge¬ 
ometrical cross-section over dust mass (see text), exemplify¬ 
ing the difference between the equilibrium radius, rgq (solid line; 
high ^), and the minimum radius, r^in (dashed/dotted line; low 
versus time t for dust approaching the star radially. Both curves 
assume = 2. Black curves represent dust entering the Hit 
region at = Myr, while grey curves correspond to the same 
situation at = 2 Myr. 


In Paper 1, it was shown that the radiation pressure force of cr 
Ori AB stratifies grains ahead of the star in a bounded region at 
0.1 pc < J < 0.4 pc for a given range of The analysis was 
limited to grains with = 1, from which a size distribution 
of the grains contained in the dust wave was inferred assuming 
compact particles, i.e., with md dust mass, a the 

grain radius and ps the specific density. However, given Eq. 
we can obtain the same values of by allowing for variation of 
2rp, and a corresponding change in For example, a grain with 
2rp = 1 and a compact geometry will have exactly the same value 
of Krp as a grain with = 0.1 and ^ increased by a factor of 10 
(a fluffy/porous geometry). We attempt to break this degeneracy 
by comparing our modeling results with observations (Sec.|^. 

The inner boundary of a dust wave is determined by grains 
that have a low ratio of Considering grains at same smaller 
values of ^ correspond to more compact grains. Due to inertia, 
these grains will slowly accelerate and, moreover, tend to over¬ 
shoot Teq due to their large intrinsic momentum, reaching the 
minimum radius rmm constrained by ^ (Fig.[^. 

For sufficient high the grains wil be momentum coupled 
with the gas of the flow. In this case, a negligible amount of 
momentum is used to accelerate the grains, from which follows 
^rad = ^drag (EQ-I^- Then, all momentum absorbed by the grains 
will be passed to me gas, and the grains will have a steady drift 
velocity throu gh the gas independent of grain size ( |Kwok|1975t 


Tielens 1983[ ). The equilibrium radius, rgq, is populated by grains 


that are momentum coupled, whose exact value will only depend 
on the flow parameters, Vf and n, and the value of the radiation 
pressure opacity 2rp- High values of ^ either corresponds to very 
small dust grains (e.g., VSGs), or large grains with a low specific 
density, such as porous grains or fluffy aggregates. It is good 
to note here that porosity and fluffiness are not directly related. 
Whereas the porosity is a measure of the volume filling factor of 
a grain (e.g., |Bazell & Dwek|1990| ) and contains no information 
on its morphological structure, a fluffy or filamentary grain is of - 
ten described by a fractal dimension (e.g., |Fogel & Leung|1998 ). 
Grains with the same porosity can have very different fractal di¬ 
mension. However in this work, only ^ enters Eq.[^ and we can 
not distinguish between porosity and fluffiness. 

Figure depicts the difference between rmin and ^eq (2rp — 
2). The curve for rmin is calculated for a 1000 A compact silicate 


grain (log(^) = 4.03). In contrast, the curve for rgq is calculated 
for a grain of similar mass, but with a geometrical cross-section 
(Td increased by a factor of 100 (log(^) = 6.03), to mimic a highly 
porous or fluffy grain that is momentum-coupled with the gas. 
The equilibrium radius rgq shrinks with t, because the density 
and velocity of the ionized gas increases, pushing the dust to¬ 
wards the star as Fdrag rises. We plot trajectories for dust ema¬ 
nating from the LI630 molecular cloud at = 0 Myr, at which 
point Fdrag is still relatively small (Fig. 0: the compact grains 
overshoot significantly, and the difference rgq - rmm is high. The 
effect of inertia is also seen in the difference in time t at which 
the grains reach rgq and rmm: the compact grains are accelerated 
more slowly, thus reaching rmm at a later stage compared to the 
grains that are coupled to the gas and reach req. However, the 
importance of the inertia of the grains will decrease with t as 
cr Ori approaches the cloud and the ionized gas flow develops 
in strength, exemplified by the curves for dust launched at to = 
2 Myr. Note that the damped oscillation of rmm around rgq is a 
peculiarity of the streamline approaching the star radially. In re¬ 
ality, dust along this streamline will flow around the star when 
the dust has reached rmm ( [Ochsendorf et al.||2014b| ) and, there¬ 
fore, only the dashed part of rmm is physically relevant. The same 
applies to momentum coupled grains: after the grains reach req, 
they will be pushed past the star. The thin part of the solid line 
reveals rgq for all momentum-coupled grains that arrive at rgq 
between 1.1 Myr < r < 2.8 Myr. 

In summary, for given flow parameters and adopting a dust 
size distribution, dust grains will be stratified in a region ahead 
of the star, stretching from req (small or porous grains) towards 
^min (compact, large grains), depending on These radii mark 
the location where the dust reaches minimal velocities with re¬ 
spect to the star, resulting in high dust densities: a dust wave. 
The radiation pressure efficiency, effectively offsets req and 
^min with respect to the star, as only enters Fj-ad in Eq. 
Values for depend on grain properties and the incident ra¬ 
diation field ( |Li & Draine||2QQH|Draine||2011| ). For a radiation 
field typical of that of an 09.5V star, compact, spherical grains 
of size ~ 1000 A can have an efficiency factor as high as ^ 
2. In contrast, small grains of 10 A have 0.1, while for the 
largest grain sizes with Ina » T, approaches unity. 


6. Modeling results 


For the flow parameters estimated in Sec. 5.2 and = 1.5, 
which is a reasonable estimate for classical silica te and graphite 
grains larger than ~ 0.01 pm (see [Draine 2011| Fig. 7), at ^ = 
3 Myr (the present situation) the equilibrium radius equals rgq 
= 0.7 pc. At the same time, rmm = 0.3 pc for 0.25 pm grains, 
typically the largest grains ava ilable in the size distribution of 
the diffuse inter stellar medium ( [Mathis et al.| 19771 [Weingartii^ 
& Draine|2001 ). Thus, for = 1.5 and assuming an ISM size 


distribution, dust will be stratified between 0.3 pc < r < 0.7 pc 
ahead of the star, which is to be compared to the observed pro¬ 
jected distances of the dust waves. Surely, the actual difference 
between cr Ori and the dust waves could be larger than the pro¬ 
jected distance due to an inclination of the system. However, this 
effect should be small, as observations indicate that cr Ori AB 
and the Horsehead (and, therefore, the LI630 molecular cloud 
associated with the Horsehead) lie approximately at the same 
distance ( [Habart et al.|2005[ ). 

In contrast to the region bounded by rmin and rgq as de¬ 
scribed above, the IR emission from the IC 434 region (Fig. 
revealed two clearly separated morphologies with corresponding 
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6.1.1. Dust population A 

The ODW at J = 1 pc places a tight constraint on the nature 
of grains responsible for morphology A. We have demonstrated 
that for the flow parameters in IC 434 (Sec. 5.2) and =1.5, 
Teq = 0.7 pc. Therefore, we conclude that in order to place dust 
at the ODW, the grains have to be momentum coupled, for these 
grains account for the outer regions of the dust wave measured 
from the star. Subsequently, we can search for the value that 
places Teq at J = 1 pc in Fig. momentum coupled grains with 
radiation pressure efficiency 2 < < 5 could reproduce the 

ODW around 1 pc. 


The required 2rp is larger than commonly accepted (Li 


|& Drain^ |2001| ). Indeed, the above described constraints 
(momentum-coupled, high Qrp), combined with the results from 
Sec. 4.3.2| (low T, insensitive to Go), seem to be incompatible 
with any gra i n material used in t he dus t models of, e.g., |Drain^ 
& Li| ( |20Q7| ); |Compiegne et al.| ( |2011| ). Here, we note that the 


Fig. 9. Location of the equilibrium radius req in pc (contours) 
of momentum coupled grains at ^ = 3 Myr (present situation), 
evaluated for a grid of radiation pressure efficiencies 2rp and 
ratios of geometrical cross-section to the mass of the grain, 
The thick black contours represent rgq = 0.1 and rgq = 1.0, re¬ 
spectively, corresponding to the peak intensities of the inner dust 
wave (IDW) and outer dust wave (ODW). The light-shaded ar¬ 
eas mark the boundaries in projected distance at which the dust 
waves are observed Fig. The black area represent grains that 
overshoot req to reach rmm, and are thus not momentum coupled. 
The vertical axis on the right denotes the grain radius if the grain 
mass is distributed in a compact geometry. 


dust arcs at a projected distance of J = 1 pc (the ODW; morphol¬ 
ogy A) and J = 0.1 pc (the IDW; morphology B). While in our 
model, Tmin can be lowered by introducing large grains (a > 0.25 
yum) into the size distribution (these have lower rmm due to their 
large momentum), req determines the outer region of the dust 
wave, which can only be increased to match the observed value 
(1 pc) by allowing for a larger value of (Eq.|^. However by 
increasing consequently, rmm will rise and divert from 0.1 
pc. Therefore, a continuous interstellar size distribution will not 
reproduce the two-folded nature of morphology A versus mor¬ 
phology B, as described in Sec.|^ and the segregation of the dust 
can only be reproduced by invoking a bimodal dust distribution, 
separating the dust into their current morphologies and observa¬ 
tional characteristics. 


6.1. A bimodal dust distribution 


grains need to be large to account for the cold temperatures (Sec. 
|4.3.2l), yet in order to assure momentum coupling with the gas, 
log(^) > 6.3 (Fig.|g. This implies that a large silicate grain with 
a mass equal to a compact grain with a = 1 jum and ps = 3.5 
g cm“^ will need an increase of geometrical cross-section by a 
factor of ~ 50 to assure momentum-coupling with the gas. At a 
= 0.1 pm, this factor decreases to 17. We defer a further discus¬ 
sion on the remarkable properties characterizing morphology A 
to Sec. Ill 


6.1.2. Dust population B 

The presence of the IDW at 0.1 < J < 0.4 pc, exhibiting the 
observational characteristics described in Sec. |4.2| indicates the 
presence of dust grains different in nature compared to morphol¬ 
ogy A. Here, there are two scenarios to consider. 

First, the Spitzer/IRS LL spectra, shown in Fig. revealed 
that the 24 pm emission from the H ii region outside of the IDW 
is dominated by carbonaceous VSGs (Sec. 4.3. 1| ). Within this 
framework, the IDW can be formed by these small, compact 
grains, that are easily momentum-coupled to the gas, yet are less 
efficient in absorbing the radiation pressure from the star. Figure 
shows that momentum-coupled grains with radiation pressure 
efficiency 0.03 < Qrp < 0.50 could reproduce the IDW at 0.1 < 


d < 0.4 pc. These values for Qrp are consistent with values ob¬ 


Section 5.4 discussed that momentum-coupled grains reach r^q, 
the outer radius of a dust wave, and that Qrp effectively offsets 
the interaction region between req and rmm from the star (Eq.[^. 
In Fig we visualize this, and plot r^q at ^ = 3 Myr for a grid 
of Qrp and This diagram reveals the necessary grain properties 
to place momentum-coupled grains at the observed positions of 
the ODW and the IDW, respectively. Above a certain value of 
the grains will obtain sufficient inertia such that they are no 
longer momentum-coupled. The division between both regimes 
shifts to large values of as Qrp increases. This is because ra¬ 
diation pressure becomes more efficient, prohibiting the grains 
from overshooting r^q. 


tained for classical compact grains of size a < 0.01 pm (Li & 
|Draine|200U |D raine 1201T1 ). 

Second, large, compact BGs that overshoot r^q due to inertia 
can explain the presence of the IDW (this scenario was already 
descri bed in Paper 1). For large, compac t grains, Qrp approaches 
unity ( [Draine & Li||2007[|Dmine||2011| ). Then, for the flow pa¬ 
rameters adopted, grains of size a = 0.5 pm reach the inner parts 
of the dust wave at J ~ 0.1 pc (in Paper 1, we de rived a = 0.33 
pm due to the different flow parameters; Sec. 5.2). This is a very 
large size above the typical cut-off for grains in the ISM (a ^ 
0.25 pm; Weingartner & Draine ( 2001| )). The exact dust grain 
size to reach the inner part of the wave depends on the adopted 
flow parameters, as well as the momentum of the grain. 

The grains that populate the IDW can therefore be explained 
by either VSGs with low Qrp that have reached rgq, or large, 
compact BGs for which Qrp reaches unity. The reduced UV-to- 
NIR opacity (Sec. 4.3.4 ) and the observe d du st color temperature 
of the IDW at 0.1 pc (73 K - 82 K, Sec. |4.2| ) hints towards large 
grains (Tab. |^, consistent with the conclusion from Paper 1: a 
carbonaceous VSG of size 1 x 10“^ pm will reach T^q = 170 K 
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at a projected distance of 0.1 pc towards the star, while a large 
grain of 0.5 yum reaches a temperature Teq = 75 K. 


7. Discussion 

The simultaneous presence of two clearly separated dust mor¬ 
phologies inside IC 434, with separate dust waves at projected 
distance of 0.1 pc and 1.0 pc, can not be explained with a con¬ 
tinuous size distribution like that observed for selected regions 
of the diffuse interstellar medium. Combining modeling (Sec.[^ 
with the constrains derived through observations (Sec. |4.2| ), we 
have argued that the dust in IC 434 is bimodal in nature. 


7.1. Derived dust properties 

Population A is characterized observationally by a cold tem¬ 
perature (20 K - 26 K; depending on the spectral emissivity J3) 
that is below the equilibr ium temperature of class ical dust grains 
predicted by DUSTEM ( [Compiegne et al.|2011| ). Moreover, the 
temperature of this population as deri ved fr om a MBB fit to the 
SED only increases by ~ 15% (Sec. |4.3!^ , while the incident 
stellar flux increases by an order of magnitude: DUSTEM pre- 
dict s a temperature increase of 50% for typical ISM grains (Sec. 
|4.2| ). We find a discrepancy between UV and IR inferred mass 
of dust contained in the ODW (Sec. |^. Modeling of the photo¬ 
evaporation flow reveals that the corresponding dust wave at 1.0 
pc can only be reproduced if the grains are momentum-coupled 
with the gas, which can be obtained by increasing the geometri¬ 
cal cross-section of the grains by factors of 17 and 50 for grains 
with a mass equivalent to a compact (silicate) grain of 0.1 jum 
and 1.0 jiim, respectively. Eurthermore, they need to be hig hly 
effective in absorbing radiation pressure (Q^p > 2; Sec. 6.1.1). In 
contrast, population B is characterized observationally by mid- 
IR at ~ 24 emission throughout the H ii region, complemented 
with the IDW that peaks in intensity at 0.1 pc. The emission in 
the Hii region outside of the IDW is caused by VSGs out-of- 
equilibrium with the radiation field, supported by the IRS spec¬ 
tra. However, the dust wave at 0.1 pc can be explained both by 
VSGs, but also by BGs that are invisible throughout the Hii 
region and only start to emit at mid-IR wavelengths close to 
the star as the grains heat up (consist ent wi th the IRS spec¬ 
trum of region 7). The modeling in Sec. 6.1.2| reveals that either 
momentum-coupled grains with Q^p < 1, or very large compact 
grains with = 1 can cause the dust wave at 0.1 pc. The color 
temperature of 73 K - 82 K, and the reduced UV-to-NIR opacity 
argue for the latter scenario (consisted with Paper 1). 

Population B is cons istent with grain properties used in cur¬ 
rent d ust models (e.g., Li & Draine 200 ![ Compiegne et al.| 
201 1| ), and we connect this population to the warm dust ob- 


modeling ( [Qssenkopf & Henning |1 994 t|Qrmel et al.|2009| ) and 
observational evidence ( [Steinacker et al.||2010| ) indicate that 
grains reach large sizes inside dense cores (> 1 yum), from which 
the presence of large aggregates seems intuitive. The observa¬ 
tions may reflect their past evolution in the cloud as the grains 
are evaporated into the H ii region (see for a similar conclusion 
Martin et~SI]|2012 ). However, the constant temperature of pop¬ 


ulation A in a varying radiation field, derived from the MBB 
fits, will pose a problem to every available dust model. A possi¬ 
ble explanation is that the size distribution is heavily affected 
inside the Hii region through processing by the intense radi¬ 
ation field of cr Ori AB. In this scenario, the radiation field 
destroys the smallest grains, which suppresses the small end 
of the size distribution. This effect was already seen in Sec. 


4.2 where a population of VSGs did not seem to survive the 


transition from the well-shielded PDR of the Horsehead neb¬ 
ula into the H ii region. An alternative explanation stems from 
the results from [ Planck Collaboration et al.| ( [^14a| ) and |Plank| 
Collaboration et al.[p014| ), who show that the peak wavelength 


of the SED does not trace the intensity of the radiation field only 
(Sec. |4.3.2| ): the authors argue that this may be a signature of 
grain properties that are evolving as the grains transit from one 
region to another. This argumentation could in principle be ap¬ 
plied to the transition of the dust in IC 434 between the dense 
molecular cloud phase to the H ii region. Still, the derived (high) 
2rp for this population remains unexplained: it is questionable 
if dust a ggregates are indeed ef ficient radiation absorbers. Eor 
example, [Bazell & Dwek| ( p^9Q| ) conclude that the optical prop¬ 
erties of fractal grains do not change significantly compa red to 


their compac t counterparts, which would imply < 2 (|Li & 


|Draine|[20Ql[ |Draine|[2()ll| ). The detailed optical properties of 


served in many H ii regi ons (|Chini et ak]! 19861 |Paladini et al.[ 
|2012t |Qchsendorf et al.||2014b| ). However, we note that Paper 
1 argued that the charging of these grains deviate significantly 
from that predicted by theory (e.g., Draine & Salpeter| 1979| ), 
consistent with the dust waves in the RCW 120 and RCW 82 
Hii regions ( [Ochsendorf et al.||2Q14b| ). However, the character¬ 
istics of population A are perhaps even more remarkable, and a 
critical review seems appropriate to place our findings in context 
with recent developments in our understanding of dust modeling 
and evolution. 

The cold temperature of morp hology A (Sec. 4.3.2| ) and the 
reduced UV-to-NIR opacity (Sec. |4.3.? ) argue for large grains, 
and the requirement for momentum-coupling with the gas in¬ 
dicates a fluffy or porous structure (Sec. |6.1|). Indeed, detailed 


grains containing chemical inhomogeneities and differing frac¬ 
tal dimensions remain poorly understood; in this respect, more 
theoretical research that constrain the optical properties of dust 
aggregates is necessary. 

In conclusion, while the analysis of the dust budget in the 
IC 434 region provides insight into the properties of the grains 
residing in the ionized gas, some of the results are incompatible 
with current modeling and should be subject to further question¬ 
ing. Indeed, evidence is accumulating that the ’standard model’ 
of interstellar dust is incapable to explain the variations of grain 
properties when moving from one region of the ISM to another 
(for a recent review, see |Jones|2014| ). The findings of this work 
support this view, and indicates that our knowledge of dust may 
be far from complete. Perhaps, we are tracing unknown grain 
materials. Eurther laboratory experiments on cosmic dust ana¬ 
logues and theoretical calculations on the optical properties of 
large aggregates must be used in designing the next generations 
of dust models, which should also incorporate the evolution of 
grain materials and properties during the lifecycle of the grains 
in the ISM. 


7.2. Limitations of the model 

The derived dust parameters depend on the adopted stellar pa¬ 
rameters, as well as the model of the approach of cr Ori towards 
the LI630 molecular cloud. Luckily, the density structure of the 
flow is well constrained by observations and reveals an expo¬ 
nential dependency with distance from the cloud. However, our 
choice fo r the ionizing flux of cr Ori AB, log(2o) = 47.56 pho¬ 
tons s“^ ( Martins et al.||200^ , is low compared to values pre¬ 
dicted by other model atmospheres. Eor example, the model at¬ 
mosphere from jSchaerer & de Koter| ( p~997| ) has log(2o) = 48.25 


photons s ^ for a 09.5V spectral type star. Because we con- 
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sider a constant effective flow thickness ( [Henney et al. 


2005), 


a higher value of Qo would only increase the density of the ion¬ 
ized gas, thereby increasing Fdrag- An increase of ionizing flux 
would therefore lead to an increase of our derived values for Q^p, 
and vice versa. 

The dependency on the flow scale height Rq is less straight- 
forwar d, an d we opted to solve this numerically. As discussed 
in Sec. |5.2[ the adopted value of Rq (= 0A6Rs) is ap propriate 
for a photo-evaporation flow emanating from a flat IF ( |Henney| 
et al.|2005] ). This choice of Rq reproduces the Her emission mea¬ 
sure with constant scale size /. In Fig we vary Rq and plot 
the dependency on our results for Qrp. Only when Rq is smaller 
than 0.46/?s, corresponding to a diverging flow, are our results 
for 2rp significantly affected. This is because in this case, the 
gas reaches higher velocities, yet has a lower density. The total 
amount of momentum transfer between the dust and gas in time 
decreases, lowering F^rag and, therefore, lowering to place 
re q at the required distances . This result was already observed 
in [Ochsendorf et ^ ( |2014b ): momentum transfer between gas 
and dust is most efficient in slow-moving, high-density photo¬ 
evaporation flows. In this respect, at high /?o, the flow parameters 
do not change significantly and barely affects Q^p anymore. 

To conclude, the observed Ha emission measure is well de¬ 
scribed by an exponential power law, n oc exp[-r//?o], which 
implies that that the gas does no t accelerate effi ciently like that 


in a globule flow, where nocr ^ (Bertoldi 1989|. Therefore, low 


values for R^ (< 0.46/?s) are unlikely for the photo-evaporation 
flow in IC 434, while our results for will not be significantly 
altered for Rq > 0.46/?s. 


7.3. PAH emission from the ionized gas 


In Sec. 4.1 we have shown that morphology A is both traced 
by thermal emission of large dust grains at T > 160 /im, and 
by IR fluorescence of PAHs, emitting at 8 jim - 12 jim. In 
the framework of dusty photo-evaporation flows, it seems pe¬ 
culiar that PAHs and large dust grains are spatially correlated. 
PAHs are easily dragged along by a flow of ionized gas, whereas 
large grains will lag behind due to inertia. We propose that PAH 
molecules are formed inside the Hii region through shattering 
upon grain-grain collisions between the separate grain popu¬ 
lations in the flow, that couple different to the gas and move 
at non-identical velocities. Sputtering is not efficient for the 
physical conditions inside the IC 434 region (Paper 1). Once 
being formed, the PAHs are destroyed quickly through photo¬ 
dissociation after UV ph oton absorptio n, expected at the inten¬ 
sity of the radiation held ( Tielens|2005| ). This could explain why 
PAHs and large grains are seen co-spatial, and results in the ab¬ 
sence of PAHs within ~ 1 pc from the star. This cavity is likely 
seen in the evolution of the PAH spectrum (Fig.|^, where a dis¬ 
continuity is observed in the 7.7 /im over 11.3 //m intensity ra¬ 
tio between region 6 and region 7. A detailed model describing 
the collision- and fragmentation process, and the evolution of 
the PAH emission spectrum inside the H ii region is beyond the 
scope of this paper, and will be presented in a forthcoming work 
(Ochsendorf et al., in prep.). 


8. Conclusions 

We have carried out a combined observational- and modeling 
study of the IC 434 region and the GS206-17-rl3 shell that sur¬ 
rounds it. We characterized the dust emission from the H ii re¬ 
gion (Secs. with the use of an updated flow model (Sec. 


Diverging flow Converging flow 



Fig. 10. The influence of the adopted structure for the cham¬ 
pagne flow on the derived values of the radiation pressure ef¬ 
ficiency, 2rp- Plotted is 2rp versus the scale height of the flow, 
Rq, normalized to the distance star - cloud, R^. The dashed line 
marks the scale height, corresponding to a flat ionization front, 
and is the one assumed in our calculations. Smaller values of Rq 
correspond to increasingly diverging flows, characterized by a 
fast acceleration of the gas. Larger values of Rq correspond to 
thick flows that accelerate more slowly. Plotted is Q^p such that 
the equilibrium radius equals rgq = 1.0 (solid lines) and rgq = 0.1 
(dashed lines). The upper curve s correspond to an ionizing flux 

, the 


(dasned lines). I lie upper curve s correspond to an ionizing ] 
of log(2o) = 48.25 photons s“^ ( jSchaerer & de Koto 1997j ), 


lower curve s correspond to an io nizing flux of log(2o) = 47.56 
photons s“^ ( Martins et al.|2005 1. 


that resulted from an analysis of the large-scale structure of 
the Belt region seen to date (Sec.[^. The important findings are 
listed below: 


1. The GS206-17-rl3 shell is driven by the Orion OBlb asso¬ 
ciation, rather than the a Ori cluster, derived through a study 
of space motions, and by comparing the kinetic energy of the 
expanding shell to the mechanical output of the Orion OBlb 
members and the weak-wind system of cr Ori AB (Sec. [^. 
The space motion of cr Ori AB is directed toward the LI630 
molecular cloud, and can be traced back to the Orion OBlc 


association, in line with the suggestion of Bally (2008) that 
assigned the cr Ori cluster to the younger Orion OBlc asso¬ 
ciation. This would solve the apparent mismatch between the 
age of cr Ori (~ 3 Myr ; |Caballero|200^ to the Orion OBlb 
association (> 5 Myr; [Bally 2008[ ). Besides that the afore¬ 
mentioned analysis of the Belt region results in a flow model 
necessary to derive the grain properties in the H ii region, it 
provides new insight in the history of the Belt region that will 
be of great interest for research on the Orion OB association 
and stellar feedback in OB associations. The Orion region 
is far from fully understood, and with the upcoming results 
from the Gaia mission, its study provides a benchmark case 
that has implications for the evolution of OB associations in 
general. 
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4. 


The global structure of the IR emission revealed two separate 
morphologies (morphology A & B) associated with a dust 
wave at 1.0 pc and 0.1 pc, respectively. We have attributed 
this segregation to a bimodal dust distribution (population A 
& B), and have characterized both populations through their 
observational properties and through modeling of the grain 
trajectories. 

A dust wave stratifies dust grains in a region stretching 
from Teq (momentum-coupled grains) towards 


(non- 

momentum-coupled grains), depending on the geometrical 
cross-section over dust mass, The radiation pressure effi¬ 
ciency, 2 rp, effectively offsets rgq and rmin from the star. 
Population A has a cold dust temperature (20 K - 27 K; de¬ 
pending on the spectral emissivity) derived from MBB fits, 
below the equilibrium temperature of ’classical’ grains at 
the incident radiation field (20 K - 55 K for grains 0.1 yum 
- 1.0 yum) predicted by DUSTEM. PAH emission is seen co- 
spatial with population A, and the increasing charge state of 
the molecules reveal that both the PAHs and population A re¬ 
side within the ionized gas. The dust wave at 1 pc, associated 
with population A, can only be reproduced if the grains are 
momentum-coupled with the gas (log(^) > 6.25), for which 
the grain need to have a significant increased geometrical 
cross section compared to compact grains. We have argued 
that the porous or fluffy aggregates that are formed inside the 
LI630 molecular cloud are good candidates to explain the 
above described characteristics of population A. However, 
the high radiation pressure efficiency (Q^p > 2 ), and the evo¬ 
lution of the temperature that does not seem to trace an in¬ 
crease of radiation field only, can not be explained by cur¬ 
rent dust models, but confirms recent results by, e.g., Planck 
([Planck Collaboration et al.|2014a| [Plank Collaboration et ^ 

|2014| ), indicatin g that our un derstanding of interstellar dust 
may be limited ( Jones|2QT^ . 

Population B is reminiscent of w arm dust, often observed to¬ 
wards the centers of H n region s ( Chini et al.|1986[ Paladini 


jet al.||2Ql^ [Ochsendorf et al.||2Q14b| ). Here, we reveal that 

it is characterized by mid-IR emission around 24 yum from 
VSGs throughout the Hit region (region 1 - 6 ), and a dust 
wave at 0.1 pc from the star where the dust reaches TE 
with Td ~ 73 K - 82 K. The dust wave at 0.1 pc can be re¬ 
produced if the grains are either momentum-coupled small 
grains (log(^) > 6.9, < 1), or large, compact grains (a ^ 

0.5 yum, 2 rp = 1) that approach the star closely due to iner¬ 
tia. These values for are consistent with prediction from 
theory. The reduced UV-to-NIR opacity, and the observed 
dust color temperature of the IDW are consistent with the 
BG scenario, confirming the results from Paper 1. 

The presence of PAH emission inside the H ii region is linked 
to shattering of grain material upon grain-grain collisions. 
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